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We study structure formation in phenomenological models in which the Friedmann equation 
receives a correction of the form i/"/r^~°, which realize an accelerated expansion without dark 
energy. In order to address structure formation in these model, we construct simple covariant 
gravitational equations which give the modified Friedmann equation with a = 2/n where n is an 
integer. For n — 2, the underlying theory is known as a 5D braneworld model (the DGP model). 
Thus the models interpolate between the DGP model (n = 2, a = 1) and the LCDM model in 
general relativity (n ^ c>o,a ^ 0). Using the covariant equations, cosmological perturbations are 
analyzed. It is shown that in order to satisfy the Bianchi identity at a perturbative level, we need to 
introduce a correction term in the effective equations. In the DGP model, E^i, comes from 5D 
gravitational fields and correct conditions on E^i, can be derived by solving the 5D perturbations. In 
the general case n > 2, we have to assume the structure of a modified theory of gravity to determine 
Ef^v We show that structure formation is different from a dark energy model in general relativity 
with identical expansion history and that quantitative features of the difference crucially depend 
on the conditions on -B^j^, that is, the structure of the underlying theory of modified gravity. This 
implies that it is essential to identify underlying theories in order to test these phenomenological 
models against observational data and, once we identify a consistent theory, structure formation 
tests become essential to distinguish modified gravity models from dark energy models in general 
relativity. 



I. INTRODUCTION 

The discovery of cosmic acceleration presents a deep puzzle for cosmology P|. A conventional way to explain 
this fact is to introduce a tiny cosmological constant or dark energy in the context of general relativity. However, 
it is also possible to think that the standard Friedmann equation which determines the expansion of the universe is 
modified. There have been many attempts to modify the Friedmann equation either empirically or based on a modified 
4-dimensional action |2| and study observational constraints coming from the expansion history of the universe. 

One example of the explicit realization of the modified Friedmann equation is provided by the Dvali-Gabadadze- 
Porrati (DGP) brane- world model 01, in which gravity leaks off the 4-dimcnsional branc into the 5-dimensional "bulk" 
Minkowski spacetime at large scales. The energy conservation equation remains the same as in general relativity, but 
the Friedman equation is modified: 

-=i^^-^p. (1.1) 

The modified Friedman equation Hl.l|l shows that at late times we have H — > Hoc ~ 1/fc- Since Hq > in order to 
achieve acceleration at late times, we require Tj, > Hq^, and this is confirmed by fitting SN observations Although 
it has been shown that the DGP model suffers from serious theoretical problems, such as the existence of a ghost in 
de Sitter solutions of Eq. and the strong coupling problem 0,0, the DGP model is the simplest covariant theory 
for modified gravity which gives accelerated expansion of the universe without dark energy. In addition, the DGP 
model allows us to determine how modified gravity affects various cosmological observations other than the modified 
expansion history of the universe. It is important to stress that in the DGP model the modification to the Friedman 
equation is derived from a covariant 5-dimensional action and junction conditions across the brane. Thus it is possible 
to derive 4-dimensional covariant effective equations which govern the dynamics of gravity on the brane. 

The expansion history of the DGP model is quite different from the LCDM model. The expansion history of the 
DGP is equivalent to that in dark energy models with an equation of state w(a) = —1/(1 + il„i{a)) where D,,„{a) is 
the density parameter for matter This connection between w and Qm can be used to impose strong constraints 
on the models. If f2m(ao) ~ 0.3, the expansion history of the DGP is equivalent to that of dark energy models with 
If > — 1. Given the fact that the present SNe data favors LCDM or even prefers w < —1, the expansion history alone 
may be enough to falsify the DGP model. Indeed, Ref. 01 showed that, combining the data from Supernova Legacy 
Survey (SNLS) and the baryon acoustic peak in the Sloan Dig ital Sky Survey data, the DGP model is not compatible 
with a spatially flat universe. This is confirmed by Ref. pil | where it is shown that the same data sets exclude the 
fiat DGP model at 3a. On the other hand, Ref. [l^ also showed that if we use the 'Gold' data for supernovae, the 
flat DGP model is still acceptable at 2a level. Together with the current status of measurements for from galaxy 
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surveys, it may be too early to reject the flat DGP model from observations. However, it is true that the flat DGP 
model is in tension with the data purely from the background dynamics. 

Then it is tempting to consider models that give different expansion histories from the DGP, like the quintessence 
models in general relativity. Inspired by the DGP model, the modified Friedmann equation of the form 

_ff^^^2_8^ (1.2) 

is proposed by Dvali and Turner ■ This is the simplest generalization of the idea of the self-accelerating universe 
where there exists a single dimensionfuU parameters r^. Thus for a given a, the model still has the same number 
of parameters as LCDM. Ref. also reported the constraint on a (—0.8 < a < 0.3 at la) from SNLS data. A 
problem of these phenomenological approaches is that it is difficult to discuss the full observational consequences. If 
we modify the Friedmann equation in the background, it is natural to think that behaviour of perturbations around 
this background is different from that in general relativity. However, without specifying self-consistent models that 
implement the modification of the Friedmann equation, it is difficult to study perturbations in a consistent way. 

The aim of this paper is moderate and we do not try to construct the underlying models. Instead, we study 
how the consistency of the gravitational equations can constrain the behaviour of perturbations around the modified 
Friedmann equation (|1.2|) . For this purpose it is desirable to derive effective covariant equations which reproduce 
the modified Friedmann equation in the cosmological background. We construct the simplest possible covariant 
gravitational equations by exten ding the covariant effective equations in the DGP model derived in Ref. |l3l | based 
on the method proposed by Ref. fl3| • A proposed covariant equation is written in the form 

Gfii, ~ K^^-)''"''H^jy — E^^, (1-3) 
where ^"^H^i, is a n-th power function of k^^G^i, — T^u and 

where k| — SttG. We assume that the conservation of the energy-momentum tensor holds. V^T^y — 0. The E'^j^ 
tensor is necessary to satisfy the Bianchi identity. 

The structure of this paper is as follows. In section II, we show that, in the cosmological background, the effective 
equations reproduce the modified Friedmann equation 1)1.1(1 with a = 2/n where n is an integer. In the background 
spacetime, it is possible to satisfy the Bianchi identity without Ef^^. In section HI, it is shown that, at a perturbative 
level, we need perturbations of E^^, to satisfy the Bianchi identity. Without an underlying model, the perturbation 
of Efj^i, cannot be determined. Hence, at this point, we have to allow ourselves an assumption about the structure of 
the theory of modified gravity. In section IV, we discuss two possibilities. One possibility is to assume that Birkoff 's 
law is respected, as is discussed in Ref. . We will show that this assumption can be translated into conditions on 
E^y. In this paper, we also consider a different possibility. Since the theory contains the DGP model as a special 
case n = 2, where the conditions on E'^;^ are known p^ . the same conditions on as the DGP may be applied for 
general n > 2. Then it turns out that the resultant theory has a very similar structure for the quasi-static sub- horizon 
perturbations to the DGP model. We comment on a potential danger of having a ghost in this case as in the DGP 
model. Then in section V we discuss how different forms of the modified theory of gravity predict different structure 
formation with the identical background expansion histories of the universe. Section VI is devoted to conclusions. 

II. COVARIANT EFFECTIVE EQUATIONS FOR MODIFIED FRIEDMANN EQUATION 

We begin with the DGP model where the covariant effective theory is well known. In the DGP model, the effective 
gravitational equations are given by ^3 

G^,, = ny)'^''^I{^.-E^„ (2.1) 

where '■^■'H^i, is 

^'^n^. = - J^Mat" + ^f^T^. + ^ [ifa.pf'^^ - [f^f] 5^. , (2.2) 
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and K^2) is given by the five-dimensional Newton constant k^^^ = SttG^. i?^;^ is a projection of the 5-dimcnsional Weyl 
tensor and it is traceless. The Bianchi identity imposes 

V^-B^, = ^^V^^'^n^,, (2.4) 

while the energy momentum tensor satisfies 

V^T^, = 0. (2.5) 
In a homogeneous and isotropic background, the energy momentum tensor is written as 

n =diag(-p,P,P,P). (2.6) 

From the definition of 2^,^, Eq. H2.3|l . we get 

n=diag(-p,P,P,P), (2.7) 



where 



Then ^^^II^j, is calculated as 



p=^H^-p, P = -^{2H + 3H^) ~ P. (2.8) 



(2)n^ - J_ ^ 0. . ) f2 9) 



" 12 V (p2 + 2Pp)5'^. 



We can show that this expression satisfies V'-^-'n^^ = 0, thus we can set £^^,y = consistently in the background. In 
fact, it is known that non-zero E^^^, corresponds to the existence of a black hole in the bulk and we can consistently 
set the black hole mass zero as a boundary condition in the bulk. The (0, 0) component of the effective equations 
gives 



B'=ri\U'--^p^^ . (2.10) 

Then we arrive at the Friedmann equation, 

- = H'- (2.11) 

where we have chosen the sign so that the solution yields a late time acceleration. In the DGP model, this corresponds 
to the choice of the embedding of the brane in 5-dimensional spacetime. 

Now we extend this covariant equation to general n. We can construct a covariant expression for '"^II^i^ that yields 



4 X 3"-i \ 





{n - l)/5" + nPp" 



(2.12) 



in the background (see Appendix). Again V^^"-'n^i, = is satisfied and we can set Ef^^, = consistently. Then the 
(0, 0) component of the effective equations gives 

^2^^2(„-l) /^2_^^\ , (2.13) 



and the Friedmann equation is obtained as 



Jf2 \ 1/n 



2(n-l) 
' r. 
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= H''fp- (2.14) 



It is useful to note that using the Friedmann equation, p and P, defined in Eq (|2.8|l are rewritten as 

P = -^(i/r,)-2(-i)/", 
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III. PERTURBATIONS 

Now we consider the perturbations. In this paper, we concentrate our attention on the scalar perturbations which 
are relevant for structure formation. We take the perturbed energy-momentum tensor as 

j.^ _[ -Sp aSq^, \ 

where 

5n] = 5n^^ ^S^Stt^. (3.2) 
Let us begin with the DGP model. The perturbed '^'ll^^ is calculated as 

<5(^)n- - 1 ( . r ~ '"'f ~ ^ (3 3) 

12 1^ -2a-V'5g" 2|(p + P)(5p + p(5P|jv - (/5 + 3P)(5^v y ' ^•^-^> 



where 



5~P ^ 5p+ \6Gl 51, = 6q,, - -^5Gl (3.4) 

6P = SP-^SG^, Sn = 6n - ^SG^. (3.5) 

and we decomposed the component of the perturbed Einstein tensor as 

SG) ^ 5G^5] + 5G^ y - ^6]6G'^\. (3.6) 

Perturbations of the divergence of '^-'II^i, are calculated as 

<5(V^(2)n^^) ^ 0^ (3 7) 

5(V''(2)n^^) ^ (P+i^ {^Sp~iaH5q + kHn)^^, (3.8) 

where the energy-momentum conservation equations 

5p + 3H{Sp + 6P) + 3{p + P)$ - a-^k^5q = 0, (3.9) 



Sq + AHSq + a" 



{p + P)^ + SP- -k^Sn 
3 



= 0. (3.10) 



are used and the Longitudinal gauge metric is adopted, 

dsj = -(1 + 2^)dt^ + a^{l + 2^)5,jdx'dx' . (3.11) 

It is clear that the Bianchi identity Eq. H2.4|l cannot be satisfied without perturbations of E^,j. We can parameterize 
the scalar perturbations of E^j^i, as an effective fluid, 

-y.i, V (3.12) 



-a bPE + Sn'^^j 

Then the Bianchi identity yields constraint equations for (5i?^jy as 

5PE + iHiSpE + 5Pe) - a-^kHqE = 0, (3.13) 
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/ 2 \ 2 H 

+ 4.H6qE + a-^ UPe - -k^SiTEj ^ --r^— {dp - ma6q + k^d^} , (3.14) 



where background equations are used to rewrite p and P. 
For general n, the covariant expression ^"^n^i/ yields 



I I —np" ^5p nap"" ^Sq^, 



4 X 3"~i \ -na 



-^p"-^6q'' n \^{n ~ l)p''~^{p + P)Sp + p"-^SP^ + fn{p,P)dTr^^ 



(3.15) 



where fn{p,P) is not determined by the requirement that the background '^"'ll^i^ is given by Eq. (|2.12|l . Again, the 
0-component of the divergence of ^"^H^i/ satisfies 



<5(v'^(")n^o) = 0, 



(3.16) 



automatically, but the i-component is non-zero. Thus we have to introduce perturbations of E^^, and the Bianchi 
identity gives constraint equations for E^,y. In this paper, we will demand that fn{p,P) is of the form fn{p,P) = 
bip^~^ + 62/5"~^P like the other components. This requirement fixes f{p, P) as 



3 5\ I 3 



Uip, P) = n + - j p^-' - "-n{n - l)p^-^P. 

We have explicitly checked this formula for ?i = 3,4 and 5 (see Appendix). 
Now we can write down the perturbed Einstein equations; 



1 



1 



na{Hr 



5Gl = -kISp H 7^—:^PE, 

na[Hrc) 



(3.17) 



(3.18) 



na{Hrc) 



SG^ = K^aSq,i 



na{Hrc) 



aSq 



E,i, 



(3.19) 



where we defined 



5G^ = kIStt 



kIStte 



2(n-l)/n 



The constraint equations for SE^^i, are given by 



6pE + AHSpE - a k Sq^ = 



(3.20) 

(3.21) 
(3.22) 



6qE + 4,HSqE + a ^ { 5Pe - ^k^SiTE 



-l§^a{Hr,){n-l)a-'F, 



(3.23) 



where 



F = 



pA 



SpE — 'SHaSqE 



k'^SiiE 



pA = Sp — SHaSq. 



(3.24) 



(3.25) 



Together with the energy- momentum conservation Eqs. (|3.9|l and H3.10|l . these equations arc the basis for the analysis 
of perturbations around the modified background Eq. 12.14f) . 
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IV. BEHAVIOUR OF PERTURBATIONS 

In this section, we study the behavior of the quasi-static perturbations on subhorizon scales. With this approxima- 
tion, the constraint equations for E^j^^, become 

Sqe = 0, 5Pe - ^k^STTE = -l^a{Hrc){n - l)a-^F. (4.1) 
The key equations obtained from the perturbed Einstein equations are 

P na(gr,) \ / 5pE_\ , . 

a? 2 \na{Hr^) -l) X na{Hr^) / ' ^ ' ' 

* + 4' = k1 -r-^ a^SnE- (4.3) 

From energy-momentum conservation, the evolution equation for matter over-density is obtained as 

A + 2HA = (4.4) 

The effective equations arc not closed and we need to assume equations of state for E^i,. Thus, at this point, we have 
to assume a structure of the theory of modified gravity. Within our approximations, we need two equations of state, 
such as 

6Pe = WiHrc)SpE, 5pE = C{Hr,)k^S7rE, (4.5) 
to close the equations. In the following, wc will consider two possibilities. 

A. Birkoff 's law 

In the background, we set E^^ = 0, thus it might be natural to consider conditions given by 

SpE = SPe = 0. (4.6) 

Wc should stress that it is impossible to take Stte = due to the constraint Eq. H4.1() . It is also important to mention 
that, in the DGP model, these conditions imply the divergent behaviour of perturbations in the bulk and so are 
unphyical. With these conditions, wc can solve the constraint equations Eq. H4.1() to get 

/ fj \ na{Hrc) (l + -1 
k^67TE^-{n~l)(-^aiHr,)] . )^ " J pA. (4.7) 

Then from Eqs. H4.2|l and H4.3|l . the solutions for metric perturbations are obtained as 



fc^ ^ K4 na{Hr, 

a? ~ Yna{Hrc) - 1 

fc2 J 1 najHr,) H a{Hr,) 

a2 M 2 na{Hrc) - 1 ^ ' [na{Hr^) - 1)2 



2* = ^rTTTTTT^pA, (4.8) 

* = ^4<!-;^ ',;t:;" . -(»-i)^ .^.,rr;" .^2 )p^■ (4-9) 



Using the background equations for a dust dominated universe. 



the above solutions can be written as 



ii^^l^i^, (4.10) 
2 na{Hr^)-V ^ ' 



= ^ff'(x)pA, (4.11) 
= -^[g'(x)+3a;g"(x)]pA, (4.12) 
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where g{x) is defined from the modified Friedmann equation Eq. H2.14|l 



ml 



Then the evolution equation for over-density is given by 

A + 2H/\ = yg'[x) + ixg"{x) 



pA 



(4.13) 



(4.14) 



These solutions are the same as the results obtained in Ref. ^3 where Birkoff's law is assumed to be respected. If 
Birkoff's law is respected, the dynamics of a spherical symmetric collapsing dust shell with radius r — R{t) can be 
derived from the Friedmann equation. The time derivative of the Friedmann equation yields 



9{x) - -xg'{x) 



which is valid also for non-flat 3-space. Then Birkoff's law implies that the dynamics of R{t) is given by 



m 
m 



g{x) - -xg'{x) 



(4.15) 



(4.16) 



where 



Klpshell 
X = — tttt; — , Pshell 



UI 



3Hi ' 



(4.17) 

(4.18) 
P 

Initially, the shell is expanding just due to the expansion of the universe, so at an initial time R{t) = a{t)Ro 
and Pshell — P- The conservation of mass in the shell means that R{t) and A are related from the condition 
R^{t) Pshell = a?Rlp as 



47ri?(t)3' 

and M is the total mass contained in the shell. The over-density is defined by 

. Pshell — P 



R{t) = a(t)i?o 1 + A(t) 



-1/3 



(4.19) 



Then the equation for R(t) can be rewritten into the equation for A. By linearizing the equation, we arrive at 
Eq. limji . 



B. DGP-like model 



In the DGP model, correct conditions for E^^, are obtained by solving the 5D perturbations and imposing the 
regularity condition in the bulk [Ifii] . The conditions are obtained as 



SPe = g^/'s, SpE = 2k'^SnE, 



(4.20) 



where the former condition comes from the fact that E^^, is a projection of 5D Weyl tensor and it is traceless. Let us 
apply these conditions for general n. Solving the constraint equations Eq. (|4.1|l . we get 



SpE = 



2 


'MHre) (l + ^^) 


- 1 


3 


na(i/r.)(l + 2(^^ 


)-'] 



pA. 



Then the solutions for metric perturbations are obtained as 

U2 ,,.2 
72 



— * 



K 



1 



(4.21) 

(4.22) 
(4.23) 
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FIG. 1: The growth history g(a) — A{a)/a is shown for n = 2 (left panel) and n — 3 (right panel). From top to bottom, the 
lines show the growth history of dark energy models with identical expansion history, the Birkoff 's law model and the DGP 
model. We set the density parameter for matter today as i}mo = 0.3. 



where 

P = l- na{Hr,) + ^^j^^ ■ (4.24) 

The evolution equation for the linear over-density is given by 

A + 2HA=!f(^l + ^yA. (4.25) 

For n = 2, these results reproduce those obtained in Ref. 0- 

Linearized perturbations are described by a scalar tensor theory with the Brans-Dickc (BD) coupling given by 

^ = ^(/3-l), (4.26) 

where the scalar field originates from the scalar polarization of the graviton. Unfortunately, this shows a potentially 
serious problem of the model, fn the Einstein frame where kinetic terms for the spin-2 graviton and the scalar mode 
are diagonalized, the scalar mode has a wrong sign for its kinetic term if 

^ + ^ = ^/3<0. (4.27) 

We can show that (3 is always negative in a phyical situation, thus this indicates that we have a ghost-like excitation. 
This can be confirmed more rigorously in the DGP model if the brane is purely de Sitter spacetime. In this case, the 
condition for /3 < implies 1/2 < Hr^ and this is precisely the condition to have a ghost derived in Ref. • We will 
come back to this problem in the conclusion. 



V. STRUCTURE FORMATION 



In general, the linear growth factor is determined by Eqs. H4.2|l . (|4.3() and H4.4|l . It is manifest that the growth factor 
crucially depends on SpE and Stte- Hence the structure of the modified gravity, which is described by the equations 
of state of E^iy, is encoded in the linear growth factor. In order to have an insight into how the different structures 
of the modified gravity change the linear growth factor, we show the linear growth factor for the Birkoff's law models 
and the DGP-like models studied in the previous section (Fig.l). We also compare them with the growth factor in 
dark energy models with identical expansion histories. As seen from Fig. 2, the growth factor is different from the 
one in dark energy models with identical expansion histories. This fact can be used to distinguish between modified 
gravity models and dark energy models , , , . 
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FIG. 2: In the left panel, the expansion histories r{z) — f^^ dzH{z)~^ are shown for LCDM, n = 4, n = 3 and n = 2 from top 
to bottom. In the right panel, the growth factor g[a) — A{a)/a is shown. The solid lines are LCDM, DGP-like models for 
n = 4, n = 3 and n = 2 from top to bottom and the dotted lines show corresponding dark energy models with the identical 
expansion histories 



The linear growth factor is the basis for tests of modified gravity against structure formation observations. But one 
also needs the metric perturbations j^. In the Birkoff's law models. 



na{Hr^) ^ _ H 



na{Hrc) — 1 



i/2 {na{Hrc) - 1)^ 



pA 



On the other hand, in the DGP-like models. 



(5.1) 
(5.2) 

(5.3) 
(5.4) 



Equations (|5.2|l and H5.4|l basically determine the integrated Sachs- Wolfe (ISW) and weak lensing effects. In the 
Birkoff's law model, this equation receives an additional correction from the modified gravity compared with general 
relativity Thus the growth factor is not enough to predict the ISW effects and weak lensing effects. Fig. 3 shows 
k^{^ — '^)/a'^KlpA. Although for n > 2, the growth factor is almost identical in both models, the ISW effects and 
weak lensing effects are different due to the time variation of the functions in Eq. 1)5.2(1 in the Birkoff's law model. This 
manifests the need for full solutions for the metric perturbations in order to predict structure formation in modified 
gravity models. 



VI. CONCLUSION 



In this paper, we studied structure formation in phenomenological models in which the modified Friedmann equaiton 
is given by Eq. ((2.14(1 . In order to study the perturbations around this modified background, we constructed simple 
covariant equations which reproduce the modified Friedmann equation in the cosmological background. By requiring 
that the Bianchi identity is satisfied, a possible form of the effective equations is restricted. At a perturbative level, 
we have to introduce an unknown term iS^i, in the effective equations to satisfy the Bianchi identity. The Bianchi 
identity imposes the constraint equations for E^^,, but they are not closed. Thus the structure of the modified theory 
of gravity is encoded in the additional conditions on which are needed to close the effective equations. 

In this paper, we analyze the quasi-static sub-horizon perturbations relevant for structure formation. We consider 
two possibilities. One is a model in which Birkoff's law is respected and the other is a model in which the perturbations 
can be described by scalar-tensor gravity as in the DGP model. These models are realized by particular choices of the 
equation of state for E^^i,. We demonstrate how the growth factor depends on the assumed structure of the modified 
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FIG. 3: - *)/a^re|pA is shown for n = 2, n = 3, n = 4 in the Birkoff's law models. - *)/a^KlpA = 1 in the 

DGP-like models and dark energy models in general relativity. 

gravity. Interestingly, in models with n > 2, the two models give almost identical growth factor. However, this does 
not imply the observational consequences are the same. The growth factor is determined by the solutions for ^P, but 
the ISW and weak lensing effects are deterixiined by the combination $ — ^. We find that the two models predict 
different behaviours of <f> — 4". This manifests the need for knowledge of metric perturbations in order to address 
structure formation in modified gravity models — the parameterization of the growth factor is not enough. 

In this paper, we only consider linear perturbations. It is crucial to study non-linear perturbations since the solar 
system constraints require that the theory must be close to general relativity at solar system scales. In the Birkoff's 
law model, it is shown that below the scale r« = {vg/ HqY^^ , the deviation frorn general relativity is small and the 
solar system constraints can be evaded This is also true in the DGP model Thus we suspect the non-linear 
recovery of general relativity is quite generic in these phenomcnological models. Detailed features of the transition 
from modified gravity to general relativity depend on the structure of the modified theory of gravity. If we can extract 
an information on the transition from linear theory to non-linear theory from observations such as weak lensing and 
cluster mass function, this also provides an interesting probe of the structure of the modified gravity. The effective 
equations are valid for non-linear physics, but we have to generalize the conditions on Efj^i, to non-linear physics. It 
would be interesting to study the non-linear dynamics such as a spherical collapse based on the effective equations 
and we will leave this issue as a future work. 

The obvious outstanding issue is to find a fundamental theory that reproduces the effective equations in this paper. 
Only the underlying theoretical model is capable of fixing unambiguously the conditions for E^^ and hence the 
growth rate . Moreover, in a self-accelerating universe in the DGP model, the additional suppression of the growth 
rate compared with a dark energy model is tightly connected with the pathology of the model. In this model, the 
scalar mode of the graviton behaves like a ghost. This is the origin of the additional suppression of the growth rate, 
because the ghost gives a repulsive force and this prevents the CDM over-density from collapsing. Our analysis may 
indicate a possibility of avoiding the ghost by changing the boundary condition for because the conclusion of 
having a theory with a negative BD parameter comes from a specific condition on Efj^^ based on regularity in the 
bulk. This modification may be achieved by an introduction of a second branc in the bulk p^ . The growth rate is 
sensitive to this modification. Hence it is crucial to develop a consistent theory for the modified gravity and, once we 
have a consistent theory, structure formation measures become essential to test the theory against observations. 
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APPENDIX A: CONSTRUCTION OF COVARIANT TENSOR ("^n^^ 

In this appendix, we construct an expressions for ("^H^^^. We first construct ^"■'H^i, as a function of T^j^ and then 
replace T^j, by T)i,y. We start from the simplest case where ^"^H^^ is a quadratic function of the energy- momentum 
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tensor T^^. In this case the answer is known from the effective equations in the DGP model, i.e. Eq. (|2.2|l . but it is 
instructive to reconstruct this expression. The general form of ^^•'ll^i, is given by 

(2)n^^ - A^T^^Tl + AaTTt + A^T^pT'^PS'^^ + A^T^Sf^^, (Al) 

where T = T^. In the homogeneous and isotropic universe, we require that '^^II^q oc and '^^■'11^,^ satisfies 

V'^(2)n^^ ^ 0. (A2) 

Then '^'ll^i/ must be of the form 



^ iB2p' + B3Pp)S\ ^ 
This requirement is sufficient to determine the coefficients Ai , A2 , Ag, and A4 up to an overall normalization; 

Al = 6A4, A2 = -2A4, A3 = -3A4. (A4) 

Then ^^^H^u is given by 



(p2 + 2Pp)(5V 

where C2 = —2A4. Interestingly, the Bianchi identity V^'^-'II^o = is automatically satisfied. 
For ri = 3, the general form of '^^'ll^i^ is given by 

(3)ni: = AiT^T^T^ + A2T^T^Ti: + A^TT^T^ + A^T^Tl^ + A^Tl^T^T^Si: + Af.TT^T^S'i + ArT^S^:. (A6) 
The same requirement as Eq. (jA3|l fixes the parameters as 

Al = 6A2 - 54A7, A3 = 9Aj, A4 = -A2 + 3Aj, A5 = -6A2 + 45A7, Ae = 2A2 - I8A7. (A7) 
Then '■^•'ll^,' is determined by 



where C = 2^42 — lAA^. Again the Bianchi identity is automatically satisfied. Now let us consider the perturbations. 
Non-trivial components are calculated as 



S^^^U\ = (6A2 - 42A7)p2%, 
S^^^nlpj^j = \{12A2-81A7)P^ + {6 A2- 36 A7)pP + 3A7p^] Sir' J, (A9) 

where S^^^ITrpj, ^ is the transverse-traceless part of the component. In the main text, we demanded that 

12^2 - 8IA7 = 0. (AlO) 

This completely fixes all components except for an overall normalization and we can verify the formula Eq H3.17|l . 

We can continue the same procedure for general n. We have checked the case for n ~ 4 and = 5. We do not 
explicitly show the results because the formula is very lengthy. For n = 4, there are 12 coefficients in the general form 
and for n = 5, there are 19 coefficients. In all cases, a requirement similar to ljA3p is sufficient to determine the form 
of (")n„„ as 











(n - l)p" + nPp" 



(All) 



The Bianchi identity V'^^"''n^o = is automatically satisfied. Then we can compute the perturbations and get 

5^"^Uirr, = fn{p,P)Sn'y (A12) 
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In this paper, wc impose the constraint that 



(A13) 



This gives n — 1 conditions. We checked that, for n — A and n — 5, these conditions can be imposed consistently and 
we get 



Uip, P) = -^4(1 V + 18pP), Up. P) = -^5(25/ + 30p3p). 



These results can be summarized as 

/„(p,p) = a 



- -n{n - l)p"-2p 



(AM) 



(A15) 
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